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ABSTRACT 

We present deep Gemini GMOS optical spectroscopy of nine luminous quasars at redshifts z ~ 
0.5, drawn from the SDSS type 2 quasar sample. Our targets were selected to have high intrinsic 
luminosities (My < —26 mag) as indicated by the [Ont] A5007 A emission-line luminosity (Irrourl)' 
Our sample has a median black hole mass of ~ 10 8 ' 8 M© inferred assuming the local Mbh-c* relation 
and a median Eddington ratio of ~ 0.7, using stellar velocity dispersions c* measured from the G 
band. We estimate the contamination of the stellar continuum from scattered quasar light based on 
the strength of broad H/3, and provide an empirical calibration of the contamination as a function of 
Lronilj the scattered light fraction is ~ 30% of L5100 for objects with -Lroin] = 10 9 5 L Q . Population 
synthesis indicates that young post-starburst populations (< 0.1 Gyr) are prevalent in luminous type 
2 quasars, in addition to a relatively old population (> 1 Gyr) which dominates the stellar mass. 
Broad emission complexes around Hen A4686 A with luminosities up to 10 8 ' 3 L© are unambiguously 
detected in three out of the nine targets, indicative of Wolf-Rayet populations. Population synthesis 
shows that ^5-Myr post-starburst populations contribute substantially to the luminosities (> 50% of 
£5100) of all three objects with Wolf-Rayet detections. We find two objects with double cores and four 
with close companions. Our results may suggest that luminous type 2 quasars trace an early stage of 
galaxy interaction, perhaps responsible for both the quasar and the starburst activity. 
Subject headings: galaxies: active — galaxies: evolution — galaxies: interactions — galaxies: nuclei 
— galaxies: starburst — galaxies: stellar content — quasars: general 



1. INTRODUCTION 

Most, if not all, bulge-domin ated galaxies harbor su- 
permassive black holes (SMB H; iKormendv fc Richstond 
119951: iMagorrian et al.lll998fl . Studying the host galax- 
ies of the most luminous quasars is essential for under- 
standing the coupled evolution of SMBHs and galax- 
ies. The linked growth has been strongly informed 
by the similar redshift evolution of th e space density 
of quasars (e.g.. iBoyle fc Terlevichl fl998: Has inger et al.1 
120051: iRichards et al.ll2006¥ and that of the star formation 
rate CSFR) (e.g.. IConnollv et alj Il997t IChapman et~aT1 
|2005| ). and the correlation between black hole mass and 
bulge stellar velocity dispersion (t he Mrh-7« relation) 
found in local inactive g alaxies (|Ferrarese fe MerrittJ 
[2000t iGebhardt et alJl2000h . 

Sim ulations of gas-rich mergers (e.g.. lDi Matteo et al.l 
2005) ascribe SMBH-galaxy co-evolution to merge- 
induced starburst and quasar activity that in turn 
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outputs energy (" feedback ") reg ulating further growth 
dSilk fc Reed Il99l iFabiad fl99l . While these models 
claim to be successful in interpreting multiple observa- 
tions (such as the quasar luminosity function and the dis- 
tribution o f Eddington ratios; e.g., iHopkins "eTail [20061 : 
IShenll200"9T ). the predicted brief phases of concurrent star- 
burst and quasar activity, in spite of being crucial for the 
growth of stellar populations and SMBHs, are rarely ob- 
served. 

There is growing evidence that low-luminosity AGN 
are not mainly induced by major mergers (e.g., 
bas ed on close neighbor or galaxy lopsidedness stud- 



ies: iFuentes- Williams fc Stockei 119881: iDe Robertis et al 
1998[ iPierce et al.1 120071: lEllison et al.l 120081: iLi et al 
20081 iReichard et alj |2009|)~ For the most luminous 



quasars (L > 10 46 erg s -1 ), on the other hand, gas- 
rich mergers have long been implicated, as supported by 
the close companions and tidal tail s seen in early HST 
studies of quasar hos t galaxies (e.g.. iBahcall et aLlll997l : 
iKirhakos e t~an il999D. but the merger scenario remains 



controversial (e.g. iFlovd et al.1120041 : (Bennert et al.ll2008l : 



lUrrutia et aTlhoOcC 

In a study of > 20, 000 type 2 AGN with L [0lu] /L Q ~ 
10 5 -10 8 ' 5 (extinction corrected) at z < 0.3 selected 
from t he Sloan Digital Sky Sur vey (SDSS; lYork et all 
I2000f) . iKauffmann et all {2003) found that the host 
galaxies of AGN with higher luminosities have younger 
mean stellar ages than a control sample of inactive galax- 
ies, while those of AGN with lower luminosities have 
stellar ages similar to those of normal early-type galax- 
ies. At higher luminosities (£[Om] > 10 8 ' 5 £©), it is un- 
clear whether luminous quasars preferentially reside in 
massive elliptical galaxies with little recent star forma- 
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tion activity (e.g-.lMcLure et al.| [l99l iNolan et~aTl l2001l: 
iDunlop et all [20031. or if there is considerable recent 
or on-going star-forming activity (e.g.. iBoroson fe Okd 
1982t iHeckman et all Il997t ICanahzofc Stockton! 120001 : 
Jahnke et al.l 12004 iLetawe et al.ll2007| ). 



Studies of the stellar populations and interaction rates 
in the most luminous quasars are significantly hampered 
by the high contrast between nuclear and stellar light. 
Techniques to circumven t this problem include t akin 
off- nucleus s pectr a Ce.g. iHutchings fe Cramptonl 1199 
Nolan et al.| 120011: Watson et alj 120081: IWolf fe Shehnsl 



20081 ) and modeling and subtract ing nuclear light in 
observed on-nucleus spectra (e.g. iMagain et al.l 120051 : 
iLetawe et all l2007t Uahnke et al.l I2007T ). Here we focus 
on luminous obscured quasars, for which the obscuring 
material acts as a natural coronagraph, allowing detailed 
study of the host galaxy morphology and stellar popula- 
tions. 

Due to the redshift evolution a nd the shape of 
the quasar luminos ity function (e.g. iBovle et "all 120001 : 
iRichards et al.ll2006l ), low-redshift quasars comparable to 
luminous high-redshift quasars are quite rar e, and in ad- 
dition the obscured ones are hard to find. iReves et alJ 
(2008) have recently published a sample of ~ 100 type 
2 quasars with Lroni]> 10 93 Lq (among a total of ~ 
900 type 2 quasars with Lp m ]> 1O 83 L ). This pa- 
per presents deep Gemini optical spectroscopy of a pi- 
lot sample of nine luminous type 2 quasars with irom] 
> 10 9 3 Lq at redshift z ~ 0.5. We determine the rel- 
ative contribution of old and young stellar populations 
to probe the starburst-quasar link. We measure stellar 
velocity dispersions, allowing us to deter mine black hole 
mass es assuming the Mbh-c» relation (|Tremaine et al.l 
2002). We combine these measurements with estimates 
of quasar intrinsic luminosities to obtain accretion rates. 

Th e strict interpret ation of AGN unification models 
(e.g., lAntonuccll Il993h states that type 1 and type 2 
quasars are identical aside from our viewing angle. Here 
we test a prediction of unification by quantifying the 
level of scattered light from the accretion disk present 
in type 2 quasars, as seen so dramatically in th e po- 
larization measurement of Zakamsk a et al. (120051) and 
the pilot imaging survey of I Zakamska et all (|2006l . see 
also ICid Fernandes et all 12004). Of course, this scat- 
tered light component is a substantial contaminant to 
our estimates of ongoing star formation, and it is thus 
critical to quantify. We also address the more funda- 
mental issue of whether or not strict unification holds; 
a substantial fraction of type 2 quasars may well suf- 
fer from galaxy-scale obscuration associated with vig- 
orous star formation, making the type 2 population as 
a wh ole more biased toward star-bursting popu l ations 
(e.g. iMartmez-Sansigre et alj 120061 : iRigbv et alj 120061 : 
lLacv et al.ll2007l ). 

The paper is structured as follows. We discuss sample 
selection in jj2j and describe our Gemini observations and 
data reduction in $3] Our data analysis method and 
results are provided in jjH We present implications and 
discussion in jj5l and summarize our main conclusions in 
Sj6] Throughout we use AB magnitudes and assume a 
cosmology with tt m = 0.3, Da = 0.7, and h = 0.7. 

2. SAMPLE SELECTION: THE MOST LUMINOUS 
SUBGROUP OF SDSS TYPE 2 QUASARS 



Our targets were drawn at the high-luminosity end of a 
parent sample of ~ 900 optically-selected type 2 quasar s 
in the SDSS ([Zakamska et alj 120031 : IReves et al.M2008l ). 
The luminosity of the [Om] A5007 emission line, L[ 0m ], 
is adopted as a proxy for the intrinsic quasar luminos- 
ity. Arising from the narrow-line region, the forbidden 
line [O in] A5007 should be much less affected by circum- 
nuclear obscuration; its luminosity is observed to be cor- 
related with the broad -band continuum l uminosity in un- 
obscured quasars (e.g. IReves et~"aTll2008l) . 

We need sufficient spectral resolution and high enough 
signal-to-noise ratio (S/N) in the continuum to obtain 
robust measurements of stellar velocity dispersions. The 
SDSS spectra have adequate spectral resolution (R ~ 
2000), but the continuum S/N for most objects at the 
high- luminosity end is inadequate, typically < 3 pixel -1 . 
Based on Lip ml measured from the SDSS spectra, we 
selected a pilot sample of nine objects with Lip ml > 
10 93 L Q (corresponding to intrinsic luminosities My < 
—26 mag; see below) for deep follow-up optical spec- 
troscopy. In addition, all of our targets were selected 
to be radio-quiet to avoid complications from radio jets 
in interpreting the results, and to be at redshifts z ~ 0.5 
for proper wavelength coverage. There is no additional 
selection criterion other than the visibility of the target 
at the observing time. In particular, we did not select 
on continuum flux, so these objects are representative of 
the parent sample at any fixed Lip ml • The targets are 
listed in Table Q] in increasing RA order. 

3. OBSERVATION AND DATA REDUCTION 
3.1. Follow-up Observation with Gemini GMOS 

Optical long-slit spectra were obtained for the nine 
targets using the Gemini Multi-Object Spectrograph 
(GMOS) on Gemini-North on 12 nights between August 
2006 and January 2007 (program ID: GN-2006B-0493). 
The seeing was variable, ranging from 0.4 to 1.1 arcsec. 
The slit-width used was 0.5 arcsec, corresponding to ~ 3 
kpc (physical) at redshifts z ~ 0.5 in the assumed cos- 
mology. The spatial sampling of the detector was 0.29 
arcsec pixel -1 . The total exposure time for each target 
ranged from 1 to 4 hours (with 30 minutes per exposure) 
depending on weather conditions, and is listed in Table 
[TJ The median S/N per rest- frame 1.0 A was 22, 32, 
and 43 over the spectral ranges 3900-3960, 4250-4320, 
and 5150-5350 A. These S/N A -1 achieved in GMOS 
spectra of our targets are at least 10 times larger than 
those of their SDSS spectra. The slit was centered on 
the quasar itself and oriented to cover as many objects 
in the field as possible to observe potential companions 
and/or extended-emission line regions. 

The R400-G5305 grating was adopted with a disper- 
sion of ~ 0.45 A pixel -1 and a spectral resolution of i£~ 
1900, spanning an observed wavelength range of 5000- 
8000 A. At the redshifts of our targets, this corresponds 
to rest-frame 3300-5300 A, a range covering the promi- 
nent signatures of both old and young stellar populations 
(see jjOJl. An AO white dwarf (G191B2B) was observed 
as a flux standard. We also observed a spectroscopic 
standard star (a K-giant) to calibrate the instrumental 
resolution (see t )4.1l for our calibration approach). 

3.2. Data Reduction 
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TABLE 1 

Type II Quasars Observed with Gemini GMOS 



Target Name 
(1) 



Redshift 


R.A. 


Decl. 


r 


^[O ill] 


Exposure 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


0.4840 


00 56 21.72 


+00 32 35.8 


20.10+0.03 


a -° -0.00 


14400 


0.5550 


01 34 16.34 


+00 14 13.6 


20.72+0.04 


q R3+0.OO 


3600 


0.4213 


01 57 16.92 


-00 53 04.8 


20.42+0.03 


M ' zo -0.00 


14400 


0.5398 


02 10 47.01 


-10 01 52.9 


20.27+0.03 


q oo+O.OO 
3,3 -0.01 


3600 


0.6261 


03 19 50.54 


-00 58 50.6 


21.32+0.07 


q 7 fi+0.00 
M '' o -0.20 


3600 


0.5560 


08 01 54.24 


+44 12 34.0 


21.13+0.05 


q 53+011 
a ' o,;i -0.00 


3600 


0.4332 


08 23 13.50 


+31 32 03.8 


20.12+0.03 


n 77+0.02 
a - ' ' -0.00 


3600 


0.5293 


09 43 11.57 


+34 56 15.9 


20.07+0.02 


y - S4 -o.oo 


5400 


0.5231 


09 50 19.91 


+05 11 40.9 


18.99+0.02 


Q 45+0.26 


3600 



SDSSJ0056+0032 
SDSSJ0134+0014 
SDSSJ0157-0053 
SDSSJ0210-1001 
SDSSJ0319-0058 
SDSSJ0801+4412 
SDSSJ0823+3231 
SDSSJ0943+3456 
SDSSJ0950+0511 



Note. — 

Col.(l): Target name in the format "SDSS Jhhmraiddmm" . 
Col. (3), (4): J2000 coordinates. 

Col. (5): r-band model magnitude and uncertainty from SDSS (uncorrected for Galactic extinction). 

Col. (6): [O ill] A5007 emission-line luminosity in the form of log(L/LQ) measured over the continuum-subtracted GMOS spectra. The 
total uncertainty is calculated as the convolution of the 1-er measurement error and the systematic uncertainty estimated using the 
difference between the GMOS and SDSS observations. A comparison between SDSS and Gemini spectroscopic measurements is presented 
in Figure [2] 

Col. (7): Total exposure time in seconds. 



The reduction of the 2D spectra was performed in 
IRAF using the Gemini GMOS package. The main steps 
to reduce each science exposure before flux calibration in- 
clude bias subtraction, flat fielding, interpolation across 
the chip gaps, cleaning for cosmic rays and bad pixels, 
wavelength calibration using arc exposures taken right 
before and after observing each target, sky subtraction, 
and extraction of the one-dimensional spectrum. The 
spectra were extracted using a 5 arcsec aperture, which 
was determined from the extension of the point spread 
function (PSF). 

The science spectra were then flux calibrated using 
the photometric standard star and corrected for atmo- 
spheric extinction using the curve appropriate for the 
Gemini observa t ory w ith IDL routines as described in 
iMatheson et al.l ([2008). The spectra were shifted to 
the heliocentric frame. Different exposures of the same 
target were calibrated separately before being co-added 
since they had different air masses and observing epochs. 
The resulting 1-D spectra of the nine targets are dis- 
played in Figure [TJ 

We check the spectrophotometry of our targets 
against their SDSS spectra, of which the spec- 
trophotometry for point sourc es is better than 5% 
(|Adelman-McCarthv etal] 120081) . The general agree- 
ment between the two observations in both continuum 
and emission-line luminosities is illustrated in Figure [2l 
While on average the two measurements agree, for sev- 
eral objects they differ from each other by more than 5%. 
We have checked that the continua of the GMOS and 
SDSS spectra of each of our targets are in reasonable 
agreement, whereas the strongest emission lines [O nr] 
AA4959,5007 A differ both in shape and in amplitude 
for several objects. Since we know the [O m] is spatially 
more extended, the line-shape differences presumably re- 
flect the different spatial coverage of the SDSS fiber (3" 
diameter) and the Gemini long-slit. Indeed, 2-D spectra 
show that the spatial extent of the strong emission lines 
is usually larger than that of the continuum. We con- 
clude that the differences between the GMOS and SDSS 



measurements, if any, are most likely due to the differ- 
ence in the long-slit and fiber coverage and slit losses due 
to seeing variations. We take these differences as the es- 
timate of systematic uncertainties. The systematic error 
added in quadrature to the measurement error is taken 
as the total uncertainty for the luminosity measurement. 

4. DATA ANALYSIS AND RESULTS 

In this section, we present our data analysis meth- 
ods and results. We measure stellar velocity dispersions 
in i )4.1i estimate scattered light from the strength of 



broad H/3 in §4.2[ analyze stellar populations in §4.3[ and 
present double cores, companions and/or extended emis- 
sion line regions covered by our long-slit observations in 

MM 



4.1. Stellar Velocity Dispersion 

Several stellar absorption features representative of rel- 
atively old populations are discernable in Figure!!] These 
include Ca K A3934 A and the G band A4304 A seen in 
most of our targets, and Mgifr A5175 A and Fe A5270 A 
apparent at least in SDSSJ0157-0053. Ca H A3968 A 
is unusable for all objects in our luminous- quasar sam- 
ple due to the overlapping strong emission lines of [Ne m] 
A3968 A and He. 

We fit the data with broadened model spectra in pixel 
space to include a power-law component to account for 
possible scattered quasar light f §4.2| . and mask spec- 
tral regions containing emission lines, detector gaps, 
an d bad pixels. We u se the direct fi t ting a lgorithm 
of iGreene fe Hoi (|2006af ) and iHo et al.l (|2009f ). with a 
model spectrum constructed by a stellar template con- 
volved with a Gaussian as an approximation of the linc- 
of-sight velocity broadening function, plus a power-law 
component. We include a power-law component here 
only to measure stellar velocity dispersions; the quan- 
tification of the non-stellar continuum is discussed later 
in §4.21 In addition, the summed components are multi- 
plied by a 3rd-order polynomial in order to model the dif- 
ference in continuum shape between data and templates, 
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Fig. 1. — Flux-calibrated Gemini GMOS spectra (smoothed by a 6-pixel boxcar for display purpose) of the nine targets in this study. 
Downward arrows indicate the expected positions of stellar absorption features Ca K A3934 A, G band A4304 A, Mgifc triplet A5175 A, 
and Fe A5270 A. Ca K and the G band are discernable in most cases. The Mgib triplet and Fe, on the other hand, show no clear detections 
except in SDSSJ0157— 0053, most likely due to the contamination by adjacent emission lines from Fe and N ions and their intrinsically 
smaller equivalent widths. In almost all targets, [Oiii] AA4959, 5007 A lines show significant asymmetry and deviation from Gaussian. 



which could result from reddenin g, template mismatch , 
and/or calibration uncertainties ((Greene fe Hoi l2006aD . 
Template stars wer e drawn from the stellar library of 
IValdes et al.l ((2004), which have a spectral coverage of 
3460-9464 A with a resolution of -1 A full width at half 
maximum (FWHM). This resolution corresponds to a ~ 
15 km s -1 , much smaller than the ~ 65 km s _1 instru- 
mental resolution of GMOS and SDSS. The difference 
between the resolution of the stellar templates and the 
instrumental resolution of GMOS has been subtracted in 
quadrature from the c* measurement. We include sev- 
eral K stars with different subclasses to account for the 
primary contribution from old populations, plus an F2 
star and a G5 star to model potential young stars. 
Since the targets are at z ~ 0.5, it was impossible to 



observe a velocity template star with the identical set-up 
as the science targets. We thus adopt the following two- 
step procedure to calibrate the instrumental resolution. 
We first obtained arc-lamp exposures before and after ob- 
serving each science target, from which we measure the 
instrumental resolution. We then use the spectroscopic 
standard star spectrum, taken with Gemini, as a sanity 
check by measuring the instrumental resolution both us- 
ing our arc-lamp procedure and also from our velocity 
templates at much higher resolution. The two measure- 
ments for the spectroscopic standard star agree, lending 
support to our overall approach of using arc exposures 
to calibrate the instrumental resolution. 

Spectral fitting was performed using each single star 
as the template at first and then with a linear combi- 
nation of different types of stars to yield best-fit param- 
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Stellar Velocity Dispersions 



TABLE 2 

Black Hole Masses, and Eddington Ratios 



Target Name 


o? 
(2) 


<x? a 
(3) 




M BH 


^BH 


^Bol 


iBol/^Edd 


(1) 


(4) 


(5) 


(6) 


(7) 


(8) 



SDSSJ0056+0032 119±16 

SDSSJ0134+0014 180±21 

SDSSJ0157-0053 212±15 

SDSSJ0210-1001 346±43 

SDSSJ0319-0058 286±15 

SDSSJ0801+4412 250±43 

SDSSJ0823+3231 348±56 

SDSSJ0943+3456 358±50 

SDSSJ0950+0511 423±61 



104±13 
157±33 

229±27 182±9 



299±35 



338±336 



363±75 



7.2±0.2 




12.8 


7.9±0.2 




13.1 


8.2±0.1 




12.7 


9.1±0.2 


9.2±0.4 


13.3 


8.8±0.1 




13.2 


8.5±0.3 


9.1±0.4 


12.9 


9.1±0.3 


9.2±0.3 


13.2 


9.1±0.2 


9.6±0.3 


13.3 


9.4±0.3 




12.9 



+0.5 
-0.5 
+0.5 
-0.5 
+0.5 
-0.5 
+0.5 
-0.5 
+0.5 
-0.6 
+0.6 
-0.5 
+0.5 
-0.5 
+0.5 
-0.5 
+0.7 
-0.5 



h2(i 
-9 
h9.8 
-3.4 
-1+2.0 
-0.65 
hl.l 
-0.39 
- + 1.9 
-0.66 

u - ''-'-0.62 

u " 5y -0.32 
, + 1.0 
-0.35 
+0.36 
-0.07 



122 

4.4j 
0.90^ 
0.50^ 
0.85_* 



0.43Z 
0.09 



Note. — 

Col. (2)— (4): Stellar velocity dispersion in units of km s — 1 (corrected for instrumental resolution; 34.lt . The total uncertainty consists 
of the measurement error and the systcmatics dominated by template mismatch estimated using templat e stars with diff e rent t ypes. 
Col. (5): Black hole mass in the form of log( -Mbh /-^0 ) inferred from er* , assuming the AfBH-c* relation of Trcmainc ct al. (20021). The 
uncertainty on AfsH listed here is p ropagated from that of er* , which does not include the intrinsic scatter in the A/sH-cr* relation 
(smaller than 0.25-0.3 dex on Afa H : ITremaine et aTll2002fl . 

Col. (6): Virial black hole mass in the form of log( A4 /Mq ). It is derived based o n the width of t he bro ad H/3 line and the intrinsic quasar 
continuum luminosity inferred from -L[oiii] (see Tablc[5j using the formula of Greene & Ho (2005). The uncertainty is dominated by 
the error of the quasar intrinsic luminosity. 

Col. (7): Quasar bolometric luminosity in the form of log(L/LQ) inferred from Liq mi . See 35.1l for details. 

Col. (8): Eddington ratio. The uncertainty is estimated with error propagation from the uncertainties of A/bh and LboI, where 

, _ 4 I CM BH m p c 
l^Edd = — • 



10.4 - 
10.2 - 
10.0 - 



9.4 9.6 9.8 10.0 10.2 10.4 
'ogfLaAoisDss 



1.0 9.2 9.4 9.6 9.8 10.0 



Fig. 2. — Continuum and emission-line luminosities measured 
from CMOS and SDSS spectroscopic observations. The SDSS ob- 
servations are made from 3-arcsec-diameter fibers, whereas the 
CMOS long-slit spectra were extracted using a 5-arcsec aper- 
ture with a slit width of 0.5 arcsec. The B-band luminosity 
Lb is calculated as th at emitted between rest-frame 3980-4920 A 
(Zakamska ct al. 2003). The emission-line [O in] A5007 luminosity 
^[O in] i s obtained from nonparametric fits using the model profile 
of [O in] lines (Reyes ct al. 2008). The symbol size corresponds to 
±5% calibration errors. The continuum and emission-line measure- 
ments from the two observations agree, given the systematics due 
to the difference between long-slit and fiber coverage and seeing 
slit losses. 

eters. There are seven free parameters in a fit with a 
single-star template (the center and width for the Gaus- 
sian, the amplitude for the power-law, and four coef- 
ficients for the 3rd-order polynomial). The power-law 
index is fixed at a\ = —1.5; we found that the disper- 
sion measurements are insensitive to the assumed values 
of a\ in the range from —1.5 to —1.0. The minimiza- 
tion of \ 2 is performed using the nonlinear Levenberg- 
Marquardt algorithm implemented by the IDL package 
"mpfit" . The uncertainties are dominated by systemat- 
ics, which include template mismatch, particularly due 
to non-solar abundances in the program galaxies and dif- 
fering intrinsic widths for different spectral types (e.g. 
iHo et al1l2009f ). For active galaxies the situation is fur- 
ther complicated by emission-line contamination (specif- 



ically from Fe and N ions in the Mgifr triplet and Fe 
regions; e.g. JGreene fe Ho1l2006af) . In order to mitigate 
these effects, we first measure the local velocity disper- 
sions over three spectral regions 3700-4020 (Ca K), 4130- 
4600 (G band), and 5080-5450 (Mg-Fe) A for every case 
in which a reasonable fit can be obtained. We define a 
fit as "acceptable" if the reduced \ 2 < 8 and the model 
traces stellar absorption features reasonably well as seen 
by eye. We found acceptable fits to the G band in all 
cases, whereas we found trustworthy fits for six of the 
nine objects around Ca K and for only one object over 
the Mg-Fe region. The Ca K and Mg-Fe regions suffer 
more severely from surrounding emission lines than the G 
band. Fitting over these two regions is more difficult also 
due to the intrinsically smaller equivalent widths (EWs) 
of Mgi6 and Fe, and the narrow cores of Ca K resulting 
from host-galaxy interstellar absorption. Despite these 
difficulties, the velocity dispersions from the Ca K and 
Mg-Fe regions agree with the G band results within 1-cr 
uncertainties in all cases. We take the G band results as 
fiducial values in the following analysis. 

Spectral fittings for the G band over the spectral range 
4130-4600 A are displayed in Figure [3] for all nine tar- 
gets, and results from the fits are given in Tabled] The 
velocity dispersions listed have been corrected for in- 
strumental resolution f §3.ip . The listed total uncertainty 
of cr* contains the measurement error and the estimated 
systematic effects, which are mostly dominated by tem- 
plate mismatch. The measured stellar velocity disper- 
sions range from ~120 to 400 km/s, with a median value 
of ^290 km/s. The inferred black hole masses and Ed- 
dington ratios are discussed in S j5.ll 

The stellar and gas velocity dispersions (estimated as 
FWHM/2.35 of [Om]) are compared in Figure 2] In our 
sample there is no correlation between them; the Spear- 
man correlation coefficient is p = —0.5 with P nu n = 0.2 
for the whole sample and p = —0.2 with P nu n = 0.7 
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Fig. 3. — Stellar velocity dispersion fits over the G band A4304 A region. In each of the nine sub-panels, the normalized spectra (un- 
smoothed) are presented on the top, along with best-fit models over-plotted as gray thick curves. Residuals from the fits are displayed on 
the bottom with null residuals marked as dashed lines. The hatched regions are those containing either strong emission lines or instrumental 
chip gaps and have been excluded from the fits. The results are listed in Table [2] 

when excluding the two double-core objects (see be- 
low). Other than the two double-core obje cts, five out of 
seven objects have cr* > FWHMjo m ]/2.35. 1 Greene et alJ 
(2009) also find no correlation in a sample of > 100 
SDSS type 2 quasars with L [0m] /L Q ~ 10 85 -10 91 . 
SDSSJ0943+3456 (the object at the upper right corner 
in Figured]) has the largest gas velocity dispersion in our 
sample (FWHM/2.35 of [Om] ~ 300 km/s). While its 
2-D spectrum shows no strong evidence for extended gas, 
it is conceivable that the gas is highly disturbed since the 
narrow lines seem to have several components at different 
velocities (Figure QJ. This might be further evidence for 
interaction in addition to its close companion (see $4jJ. 

Two caveats must be mentioned here. First, roughly 
half of elliptical galaxies in SDSS with the apparently 
largest velocity disp ersions are in fact close doubles 
(|Bernardi et ai.ll2006[ ). Similarly, our measurements may 
be biased by additional stellar components covered by 
the slit. Indeed, as discussed in Sj4~I SDSSJ0950+0511 
and SDSSJ0823+3231 are both such double-core sys- 
tems marginally resolved in our long-slit observations. 
As discussed in §5.1[ their black hole masses (Edding- 
ton ratios) listed in Table [2] are therefore likely to be 
over-estimated (under-estimated). Second, our results 
have not been corrected for aperture effects. At red- 
shifts of 0.5, the size of our aperture (corresponding to 
~ 30 kpc) is large enough that the old stellar popula- 




100 300 
FWHM t0 m/2.35 (km s"') 

Fig. 4. — Comparison of stellar and gas velocity dispersions. 
The stellar velocity dispersion <r» is from the G band (Figure [3}, 
the error of which is dominated by template mismatch. The gas 
velocity dispersion is determined as FWHM/2.35 of [Om] 5007, 
which exhibits asymmetry and/or non Gaussianity for most ob- 
jects in our sample (Figure [TJ. The uncertainty of FWHMjom] 
is estimated using the difference between SDSS and GMOS mea- 
surements (also see Figure The two objects marked with open 
circles have resolved double nuclei in the projected central 5 arcsec, 
s o the measured c* might over-represent that of each component 
( 34.41 1. There is no correlation between stellar and gas velocity dis- 
persions in ou r sam ple whe ther or not the two double-core objects 
are included ( 94.11 also see Greene et al. 2009). 
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Fig. 5. — Quantifying the broad H/3 component to con- 
strain scattered light contamination of the stellar continuum. The 
two examples represent null (SDSSJ0056+0032) and significant 
(SDSSJ0823+3231) detections, respectively. Each of the narrow 
emission lines, H/3 and [Om] AA4959, 5007, is modeled using a 
template combining t he blue wing of [ O m] A4959 and the red 
wing of [Om] A5007 (Reye s et al.ll2008T) . The broad H/3 compo- 
nent is modeled as a Gaussian. Top row: spectra are presented 
over the fitting range 4780—5100 A. Data are plotted in black 
and best-fit models are in green. Gray curves show the contin- 
uum model which has been subtracted before line fitting. Middle 
row: continuum-subtracted spectra zoomed in at the base. The 
model for the broad/narrow H/3 is shown in red/blue. Bottom 
row: residuals from the fittings (smoothed with 6-pixel boxcar) 
are shown (solid), along with l-o~ error spectra (dotted) and null 
residuals (dashed). 

tions from which we measure cr* also contain those in 
the disk component, if any, so that the c* (and by exten- 
sion black hole masses) are likely to be over-estimated. 
However, we estimate that this contamination should be 
small, considering that most of the host galaxies of lu- 
minous type 2 quasars in our parent sample are likely 
to be dominated by ellipticals with de Vaucouleurs light 
profiles as seen in the imaging study of a pilot sample 
(Zakamska ct al. 20061). For those with disk components, 
the contamination should still be moderate since the old 
stellar populations are predominantly in stellar bulges. 

4.2. Quantifying Scattered Quasar Light from Broad H[3 
in Direct Spectra 

Using spectropolarimetry, I Zakamska et al.1 (|2005l ) de- 
tected high polarization levels in a sample of SDSS type 
2 quasars as luminous as our targets. Significant scat- 
tered light has also been revealed by HST imaging in 
a so mewhat less luminous SDSS type 2 quasar sam- 
ple ([Zakamska et al.l 120061 ). Early studies attributed 
the blue light frequently observed in the host galax- 
ies of obscured AGN to scat tered light predi cted by 
AGN unification models (e.g.. lAntonuccil Il993f ). How- 
ever, it was later realized that in many cases scat- 
tered ligh t explained only part of the blue contin - 
uum (e.g., lKavl[l99l iCid Fernandes fc Terlevichl[l995l) . 



as the polarization fraction observed in the contin- 
uum is smaller than that in the permitted emission 
lines, requiring an additional source of non-po l arized 
blue light (e.g.. iGoodrich fe Millerl [1981 iTranl fl995l : 
iHeckman et al.lll995ft .~ Since there is strong evidence for 
the presence of both blue starlight from young stellar 
populations and scattered quasar light in the continuum, 
we need to determine and remove the contamination from 
scattered light before performing stellar population mod- 
eling to the continuum. 

To quantify the level of scattered quasar light, we 
model the H/3 line as the sum of a narrow and an un- 
derlying broad component from scattered light. Then 
we make use of the tight correlation between the H/3 and 
the c ontinuum luminosit ies observed in unobscured AGN 
(e.g. JGreene fc Holl2005l ) to infer the corresponding non- 
stellar cont inuum. The rms sc atter of this correlation is 
- 0.2 dex ((G reene fc Ho 2005|) and is added in quadra- 
ture with the measurement error into the uncertainty of 
our estimate of the scattered light. The assumption of 
our approach is that if broad lines are present with a cer- 
tain flux, the quasar continuum should also be present at 
a level proportional to the broad emission line flux. H/3 
is adopted because it is the strongest permitted line in 
the wavelength range covered by the GMOS spectra. 

The procedure to quantify the broad H/3 component is 
the following. First a continuum model is constructed by 
performing a x 2 fit to the emission-line-masked spectrum 
over the range of 3600-5600 A using multiple in s tanta - 
neous starburst templates of iBruzual fc CharlotJ (|2003l ) 
broadened with the measured stellar velocity disper- 
sion and a power-law component assuming a\ = —1.5. 
Next we make a four-component fit (broad H/3, narrow 
H/3, [Om] A4959, and [Om] A5007) to the continuum- 
subtracted GMOS spectrum over the range of 4780- 
5100 A. We fit each of the three narrow emission 
lines non-parametrically with an [O m]-profile model con- 
structed using the blue wing of JO ml 4959 and the red 
wing of [Om] 5007 (|Reves et al.ll2008l) . For objects like 
SDSSJ0943+3456, which appears to have several compo- 
nents with different velocities in its narrow lines, we do 
not need to use multiple components to fit the narrow- 
line profile, since they are accounted for all at once using 
the model constructed from the [O m] lines. Unlike para- 
metric multi-component fits, this non-parametric ap- 
proach is free from degeneracies between multiple nar- 
row and broad components, introducing no additional 
uncertainty to the measurement of broad H/3. The as- 
sumption here is that the narrow H/3 component has the 
same profile as that of [O m] . Their rcdshifts are con- 
strained to be the same as they all originate from the 
narrow-line region. The broad H/3 is fit as a single Gaus- 
sian. There are seven free parameters in one fitting (a 
rcdshift for the three narrow lines, an amplitude for each 
narrow line, and three parameters for the Gaussian). Al- 
though in principle one should iterate the continuum and 
emission-line fits, in practice we have found the broad H/3 
component to be insensitive to the details of the contin- 
uum fit so that one it eration is enough. 

As pointed out by iReves et al.1 (|2008f ). it is impor- 
tant to use the [Om] line profiles instead of Gaussians 
or Lorentzians to fit the strong narrow emission lines. 
Almost all the objects in our luminous-quasar sam- 
ple have narrow line profiles which deviate significantly 
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TABLE 3 

Scattered Light Quantification and Continuum 
Decomposition Based on Broad-H/3 Measurement 



-"5100 ^510( 

Target Name L^ r 5 ad (10 43 erg/s) (%) (%) 
(1) (2) (3) (4) (5) 



w/o broad H/3 



SDSSJ0056 . . . 


<7.29 


<0.74 


< 9 


<0.2 


SDSSJ0134. .. 


<7.66 


<1.6 


<17 


<0.2 


SDSSJ0157. .. 


<7.04 


<0.45 


< 8 


<0.2 


SDSSJ0319 . . . 


<7.66 


<1.6 


<17 


<0.2 


SDSSJ0950 . . . 


<7.69 


<1.7 


< 5 


<0.4 


with broad H/3 


SDSSJ0210 . . . 


7 8zl+ - 21 


n Q+1.2 
z -°-1.3 


37±22 


n 2 +0 - 4 
u - z -o.i 


SDSSJ0801 . . . 


7 7 r,+0.26 

'•' u -0.79 


1 7+ 1 ' 2 


45±36 


n o+0.9 
u - a -0.3 


SDSSJ0823 . . . 


8.1b_ n 


4 Q+ - 8 
*- d -0.8 


55±12 


n 4 +0S 

u - 4 -0.3 


SDSSJ0943 . . . 




o 7 +0.9 
1.0 


28± 8 


n o+0.6 

u.a_ 2 


Note. — 



Col.(l): Abbreviated target name. 

Col. (2) : Luminosity of the broad-H/3 component in the form of 
log(L/ Lq) determined from emission-line fits over the continuum- 
subtracted GMOS spectra. For objects with null broad-H/3 detection, 
1-a upper limits arc given. The fitting method and error estimation are 
described in S|4.2| 

Col. (3): Quasar (scattered light) monochromatic luminosity L5100 — 
\L\ at rest-frame A — 5100 A, determine d from the broad-H/ 3 lumi- 
nosity, using the Lsiqq-Lhq calibration of IGreene & c Ho (200|]). The 
uncertainty contains both that propagated from the broad-H/3 luminos- 
ity and the scatter in the observed Lqiqq— Lh/3 correlation. 
Col. (4): Percentage of scattered quasar light relative to the total ob- 
served luminosity at rest-frame A — 5100 A. The uncertainty is prop- 
agated from the errors of £5100 ere an d -^5100- Since this ratio is a 
relative quantity, it is independent of the calibration uncertainty of the 
GMOS spectra. 

Col. (5): Percentage of the nuclear quasar light scattered into our line of 
sight. The intrinsic £5100 (u nobscured) is infer red from £[onr] using 
the L [o nil -M2500 calibrati on jRe yes et al. 2008) and assuming a spec- 
tral index of a v — -0.44 (Vandcn Berk et al-IPOOlJ) . The uncertainty 
contains both that pro pagated from ^^ c jnn tcrcd and the scatter in the 
^[O ill] -^2500 relation (Reyes et al. 2008). 

from Gaussian and Lorentzian profiles, often with strong 
asymmetries and sometimes double peaks. The specific 
profiles of the strong narrow emission lines, if not prop- 
erly accounted for, can cause false broad-H/3 detections. 

The broad-H/3 measurements are summarized in Table 
[3l We unambiguously detect broad H/3 components in 
four objects; we list 1-a upper limits on broad H/3 for 
the other five objects. The uncertainty on the broad H/3 
luminosity is estimated using the 1-a flux-density error 
spectrum integrated over 4794-4932 A, which is equiv- 
alent to assuming a Gaussian with a FWHM of 8000 
km/s centered at 4863 A. This is a rather conserva- 
tive estimate, considering that the broad-H/3 components 
detected in our sample all have FWHM > 4000 km/s. 
Figure [5] presents examples of significant and null broad- 
H/3 detections, respectively. The unambiguous broad-H/3 
detections or robust upper limits rely on the high S/N 
achieved by our GMOS observations (a median S/N of 56 
per rest- frame 1.0 A over the spectral range of 4841-4881 

A). _ 

Using the correlation between H/3 luminosity and 
monochromatic continuu m luminosity £5100 d etermined 
from unobscured AGN bv lGreene we quan- 



tify the amplitude of the scattered-light component 
^lioo tcred ( or ^ s upper limit when un-detected) as listed 
in Table [31 We adopt a scatt ered light spectrum with a 
spectral index of a v = -0.44 (fVanden Berk et al.ll2001i) . 

One caveat that must be mentioned here is that we are 
measuring the scattered light inside our slit, the contri- 
bution of which could vary with slit posi tion because the 
scatt ering can be very asymmetric (e.g..lZakamsk a"et al.l 
2006). In particular, if the slit is oriented perpendic- 
ular to the axis of the scattering cone, the scattered 
light covered by long-slit spectroscopy could be much 
smaller than that detected by broad-band polarimetry. 
In our sample, only SDSSJ0056+0032 has broad-band 
(observed 480-600 nm) polariz ation measurement (P . 
Smith, private communication; Zakams ka et al.l [2006) , 
which indicates a high level of polarization (10.2 ± 1.6%) 
suggesting a significant scattered light contribution (> 
10%) around rest-frame ~ 3600 A. A scattered light 
fraction of 10% at ~ 3600 A will translate into 5% 
at 5100 A (both in rest-frame) given the spectrum of 
SDSSJ0056+0032, assuming a v = -0.44, which is not in 
conflict with our 1-a upper limit (Lf 1 ^ cred / < 9%). 
However, the estimate of scattered light fraction at 5100 
A based on polarization measurement can be larger than 
5% assuming a more typical intrinsic continuum polar- 
ization (e.g., 20%), and thereby differ from our estimate 
based on the strength of broad H/3. This possible dis- 
crepancy may be due to the difference in coverage of 
the long-slit and broad-band observations. Indeed, the 
GMOS slit was oriented within 15% of the polarization 
position angle and was therefore almost perpendicular to 
the scattering direction. 

4.3. Stellar Populations 

Our Gemini spectra cover a spectral range containing 
various signposts for both young and old populations. 
These include: (1) the 4000- A break and stellar absorp- 
tion lines (such as Ca II K&H at 3934, 3968 A, the G 
band at 4304 A, the Mg 16 triplet at 5175 A, and Fe 5270 
A) indicative of relatively old populations (i.e., with ages 
> a few Gyr); (2) the Balmer continuum limit at 3646 A 
and the Balmer absorption-line series representing post- 
starburst populations (with ages < 1 Gyr); and (3) the 
broad Wolf-Rayet emission complexes around ~ 4660 A 
which trace very recent starburst activity in the last 5 
Myr. 

As foreshadowed by iHeckman et al.l (jl997) in their 
study of Mrk 477, one of the most luminous local type 
2 quasars, detecting post-starburst signatures through 
Balmer absorption becomes more difficult at high lu- 
minosity, because of the strong Balmer emission lines 
excited by the AGN. Indeed as seen in Figure [TJ the 
Balmer absorption series is barely discernable for most 
of the nine luminous quasars. The only exception is 
SDSSJ0056+0032 (which lies at the low-luminosity end 
of our sample, and has prominent very young stellar pop- 
ulations, see below) , for which the high-order Balmer ab- 
sorption lines (H8, H9, and H10) are visible. 

The broad emission complex around He n A4686 A due 
to Wolf-Rayet (WR) stars, on the other hand, is more 
prominent in luminous type 2 quasars. WR stars are 
very massive evolved stars that ap pear only ^2-5 Myr 
after a burst of star formation (e.g. lVacca fe Contilll992l ; 



HOST GALAXIES OF LUMINOUS TYPE 2 QUASARS 



9 



TABLE 4 

Measurements of Wolf-Rayet Features around He ii4686 



Target Name(l) 


SDSSJ0056 


SDSSJ0134 


SDSSJ0950 


£nv4610 


•(2) 








EWnv4610 


.(3) 






0.71 ± 0.12 


FWHM NV 4610 


■(4) 






1520 ± 280 


III 4640 


.(5) 


7.04+°-°* 

— U.Uo 


7 on+0.07 
— U.Uo 




EWn III 4640 • ■ ■ 


.(6) 


2.4 ±0.2 


3.9 ±0.7 




FWHM NIII4 640 


•(7) 


1920 ± 200 


2680 ± 540 




^Hc II 4686 


.(8) 


7 ,0+0.05 
'•18-0.06 


7 40+ - 04 
' - 4u -0.05 


o no + 0.01 

»-^»_0.01 


EWhoII4686 •■• 


.(9) 


3.3 ±0.4 


4.5 ±0.5 


11.9 ±0.3 


FWHM He ii4686 


(10) 


1830 ± 250 


1600 ± 230 


6410 ± 250 


%E 


(11) 


7.42 ± 0.05 


7.67 ±0.06 


8.31 ±0.01 


Type 


(12) 


WN6-8 


WN8 


WN3-4 


^WR 


(13) 


10 5 


10 5 


10 6 


L H{3 


(14) 


0.10 ± 0.01 


0.14 ±0.02 


0.80 ±0.03 



Note. — 
Row(l): Abbreviated target name. 

Row(2) , (5) , (8) : Luminosity of the WR feature in the form of 
log(L/L ). 

Row(3),(6),(9): Equivalent width in A of the WR feature. 
Row(4),(7),(10): Full width at half maximum in km/s of the WR fea- 
ture. 

Row(ll): Total luminosity combining all the WR features detected 
around He II 4686 in the form of log(L/_LQ). 

Row(12): Subclass of WR stars estimated based on the intensity ra- 
tio between differe nt WR features according to the classification of 
ISmith et atl (Il996lj ■ 

Row(13): Number of WR stars present estimated from the luminosity 
and the subclass. 

Row(14): Ratio of the WR and H/3 luminosities. This can be esti- 
mat ed as the low er limit to JV^r/^o to first-order approximation 
(e.g. [Srnith|[l99iT) . 



iSchaerer fc Vaccal Il998f ). whose several km s _1 winds 
produce the broad emission features. Because of their 
short duration, WR stars are an excellent clock for very 
recent starburst activity. In addition, the subtypes of 
WR stars (e.g. WN or WC) can be inferred from the rel- 
ative strengths of different WR features, and the WR lu- 
minosity can then be furt her translated into the number 
of W R stars present (e.g. JSmi"thlll991t ISchaerer fc Vaccal 
[1991) . 

4.3.1. Wolf-Rayet Populations 

In this section, we address the frequency of WR pop- 
ulations in luminous type 2 quasars. Conspicuous WR 
signatures in the optical are less affected by dust obscu- 
ration than are UV tracers in chara cterizing starburst 
activity (e.g. iKunth fc Contini| [l999). WR populations 
have been observed in a number of quasar hosts in the 
nearby universe, including the f amous infrared-luminous 
galaxy IRAS 9104 + 4109 (jKleinmann et all 119881 : 
iTran et aT1l2000D. and three luminous type 2 S eyfert nu- 
clei studi ed by iGonzalez Delgado et alj (|200lD including 
Mrk 477 (|Heckman et al.lll997l ). However, the frequency 
of their occurrence and their relation to the nuclear ac- 
tivity in luminous quasars remain important open issues. 

In order to quantify WR populations, we fit the broad 
WR emission features and the surrounding nebular emis- 
sion lines simultaneously to the continuum-subtracted 
CMOS spectra over 4550-4800 A. Here we use "neb- 
ular" to denote the emission lines that do not arise from 
WR stellar outflows (which are narrow for the forbidden 
lines and could have broad bases due to scattered light 
for the permitted lines). Fitting bot h simultaneously is 
required to isolate the WR features (Brinchmann et al. 



2008), but is only possible with high S/N spectra. The 
continuum model is constructed in an identical way as 
that in the broad H/3 fits fi P~2"|) . 

The WR features to be detected include: N v A4610, 
N in A4640, C in/iv A4650, and Hen A4686 A. They are 
fit with four Gaussians. While there are no constraints 
applied on the widths of these four WR features, only 
those broader than the narrow emission lines are consid- 
ered as detections. The WR features are broadened as 
they arise from stellar outflows and the widths can vary 
significantly among individual WR stars. 

The surrounding nebular emission lines include [Fe m] 
A4658, [Fe m] A4669, Hen A4686, [Fe m] A4701, [Ar iv] 
A4711, He i A4713, [Ne iv] A4714, [Ne IV ] A4725, [Ar 
ivl A4740, and [Fe ml A4755 (for the line lis t, see e.g., 
iHeckman et~aLlll997tlBrinchmann et al.| [2008U The neb- 
ular He ii A4686 is modeled with the profile of H/3, which 
conveniently accounts for the scattered light, if any, to 
avoid false WR detection of Hen; each of the other nebu- 
lar emission lines is fit with a single Gaussian, the width 
of which is fixed to be that from a combined fit of the 
strong narrow lines (narrow H/3 and [O in] AA4959,5007). 
Gaussian profiles are adequate for these relatively weak 
forbidden lines (from ions of Fe, Ar, and Ne, with Ar 
lines being the strongest among them), as we find that 
given the S/N of our spectra, the asymmetries in these 
lines are insignificant. The redshifts of all the features 
are fixed to be the same as that of the strong [Om] lines. 

WR features have been unambiguously detected in 
three out of the nine targets, as shown in Figure El Also 
displayed in comparison are two objects with null de- 
tections and one object with only suggestive WR fea- 
tures, which we do not count as a detection (see figure 
caption for details). We list results for the three tar- 
gets with definitive WR detections in Table |U Lumi- 
nosities emitted in each of the WR features and the to- 
tal WR luminosity are presented with measurement un- 
certainties. The WR luminosities are comparable to or 
higher than the most luminous WR galaxies k nown (e.g. 
lOsterbrock fc Cohenlll98l lArmus et al.lll988h . The in- 
ferred subtypes and the numbers of WR stars are also 
listed, along with the ratio between WR and H/3 (nar- 
row) luminosities which can be viewed to first-order ap- 
proximation as the ratio between WR and O stars (e.g., 
lSmith|[l99ll ). The three targets wit h unambiguous WR 
populations all have WN subtypes (ISmith et al.|[T996h . 
similar to Mrk 477 (jHeckman et al.lll997f ). The C m/iv 
bump at 4650 A indicative of WC subtypes is not de- 
tected in any of the nine objects. We list the estimates 
of WN subtypes from the inte nsity ratio between differ- 
ent WR features according to ISmith et al.l (l996). The 
line widths are consistent wit h those expected from the 
corresponding WN subtypes (|Smith et al.| [T996). 

4.3.2. Population Modeling of the Stellar Continuum 

We now fit stellar population models to the stellar con- 
tinuum to get a more complete view of the mix of ages in 
the host galaxies. First we fit the stellar continuum us- 
ing a linear combination of nine instantaneous starburst 
models with ages of 0.00 5, 0.025, 0.10, 0.29 0.64, 0.90, 
1.4, 2.5, and 5.0 Gyr from lBruzual fc Chariot! (|2003l ). af- 
ter masking emission lines and detector chip gaps. We 
use those nine age grids to sample the whole history of 
a z ~ 0.5 galaxy, having in mind a model in which the 
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Fig. 6. — Broad Wolf-Rayet (WR) emission complexes around Hen A4686 A. They include: N v doublet A4610 A, N in A4640A, C nr/iv 
A4650A, and Hen A4686 A. In each example shown, data (black) and best-fit model (dark grey) are displayed on the top and residuals from 
the fits (smoothed by a 6-pixel boxcar) are shown on the bottom along with l-er error spectra (dotted) and null residual (dashed). Models 
for WR features are shown in red. Plotted in blue is the model for He n A4686, the strongest nebula r line in this range, which is fit non- 
parametrically with the observed H/3-line profile. Details of our fitting method are presented in 34.3.1I Two examples of null WR detection 
are presented in the top row: SDSSJOf 57— 0053 has no detectable scattered light, while SDSSJ0823+3231 has scattered light detected in a 
broad component to H/3 (Figure^. WR features are detected unambiguously in three targets (SDSSJ0056+0032, SDSSJ0I34+001 4, and 
SDSSJ0950+05fI). The suggestive WR feature is less convincing in SDSSJ080f +4412: the apparent WR component of Hen 4686 could in 
fact be due to the degeneracy with the nebular component, and the detections of the other three WR features are marginal. Results from 
the three clear detections are summarized in Table [4] 



bulk stellar population in a galaxy was built up in mul- 
tiple starburst events at different epochs. For objects 
that have broad H/3 detections ffc |4.2p . a model of scat- 
tered quasar light is subt racted assuming a v — —0.44 
(jVanden Berk et "all l200lh . In producing the instanta- 
neous starburst model, we assume "Padoval994" stel- 
lar evolutionary tracks (|Bruzual fc Chariot 2003 ) and 
a Chabrier initial mass function (IMF; |Chabriei] |2003) 
with a low and high mass cutoff of 0.1 and 100 M@, 
respectively. We adopt solar metallicity in the baseline 
models; assuming super-solar metallicity would result in 
smaller age estimates whereas assuming subsolar metal- 



licity would result in larger ages. Template spectra are 
broadened by the measured stellar velocity dispersions 

(S3). 

Fitting was performed over the rest-frame spectral 
range of 3650-5600 A. Fitting instead over the range 
3600-5600 A results in a slightly higher contribution 
from young stellar populations for the objects with the 
largest L[oin] hi our sample. This is most likely due 
to the Balmer continuum emission blue-ward of 3646 A 
from the quasar which could bias stellar ages towards 
younger values. 

We have also tested fitting the stellar continuum us- 
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Fig. 7. — Reduced x 2 from stellar population fitting using a 
single-age instantaneous starburst model. Contamination from 
scattered quasar light has been excluded before the fitting. Re- 
sults for d ifferent targets are color-coded as labeled on the plot. 
See 34.3.2l for more discussion. 



TABLE 5 

Mass and Age Estimates of Stellar Populations 







ti 


/ £5100 


/m» 


ti 


SFR 


Target Name 


M» 


(Gyr) 


(%) 


(%) 


(Gyr) 


(Moyr- 1 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


w/o broad H/3 


SDSSJ0056 . . . 


10.0 


0.005 


53 


2 


1.1 


36 


SDSSJ0134 . . . 


10.3 


0.005 


46 


1 


2.5 


42 


SDSSJ0157. .. 


10.3 


1.4 


14 


33 


2.5 




SDSSJ0319 . . . 


10.8 


0.025 


19 


1 


5.0 


20 


SDSSJ0950 . . . 


10.9 


0.005 


64 


1 


5.0 


150 


with broad H/3 


SDSSJ0210 . . . 


8.3 


0.005 


100 


100 


5.0 


44 


SDSSJ0801 . . . 


10.0 


0.005 


38 


2 


5.0 


14 


SDSSJ0823 . . . 


10.3 


0.025 


23 


1 


5.0 


8 


SDSSJ0943 . . . 


10.2 


0.025 


80 


11 


5.0 


(58 



Note. — 
Col.(l): Abbreviated target name. 

Col. (2): Estimate of total stellar mass in the form of log(M*/AfQ). 

Col. (3): Age estimate of the younger one of the two stellar populations 

which have the highest contribution to Lsioo- 

Col. (4): Percentage contribution of the population in Column 3 to L5100 
after scattered light is subtracted. 

Col. (5): Percentage contribution of the population in Column 3 to stellar 
mass. 

Col. (6): Age estimate of the older one of the two stellar populations which 
have the highest contribution to Lsioo- 

Col. (7): Average star formation rate estimate in the past < 0.1 Gyr. 

ing only one instantaneous starburst template to the 
scattered-light subtracted spectra (using scattered-light 
models constructed with either detected broad-H/3 com- 
ponents or upper limits). The reduced \ 2 as a function 
of age is shown in Figure [7] For each object there are 
two ages at which the reduced x 2 reaches a local min- 
imum. In almost all cases, they represent two popula- 
tions, one relatively young (<0.1 Gyr) and one relatively 
old (> 1 Gyr), which dominate the starlight. Fitting a 
single-age model prefers a young population (<0.1 Gyr) 
over an older (>1 Gyr) one in most cases. Thus, the 



data cannot be explained by pure old populations plus 
scattered light; there is a substantial contribution from 
young populations. The ages of the young populations 
that we find are significantly smaller than the typical 
post-starburst ages (> 1-2 Gyr) estim ated in host galax- 
ies of quasars with lower luminosities (iKauffmann et alJ 
l2003HCanaiizo et al.ll2007t iBennert et al.ll2008t see 

Similarly, in the nine-component fittings, we find that 
two components dominate the fits in almost every case, 
one relatively old (> 1 Gyr) and one relatively young 
(< 0.1 Gyr), so that we redo the fits with just these two. 
The two components with the highest weights in a nine- 
component fit also produce the best fit when we allow 
the ages to vary (among the adopted nine age grids) in 
a two-component fit, which justifies our choice. Two- 
population fitting produces a reduced \ 2 at least 10% 
smaller than that when we fit the data with a single 
population, while adding a third population does not sig- 
nificantly improve the fits. 

Figure [8] displays two examples of population synthe- 
sis results for objects without detectable scattered quasar 
light, and the fit for the object with the highest scattered 
light fraction is illustrated in Figure [9l In Table [5] we list 
the estimated total stellar masses, ages, and fractional 
contributions to stellar luminosity and mass of the two 
populations. The uncertainty of the age estimates can 
be inferred from the population model grids we use. For 
example, if the listed best-fit age is 0.005 Gyr, the next 
model examined was 0.025 Gyr, so its uncertainty can 
be estimated as < 0.025 Gyr. We find that eight of our 
nine targets contain a very young (< 0.1 Gyr) and an old 
(> 1 Gyr) population; in half of these, the age estimates 
for the young components are < 25 Myr. The one object 
which does not have a considerable contribution from a 
< 0.1 Gyr population (SDSSJ0157-0053) has an age es- 
timate for its younger population of 1.4 Gyr. The young 
stellar populations in the three objects with Wolf-Rayet 
populations (JOT} all have best-fit ages of 5 Myr for 
the young population, and this population contributes 
> 50% of £5100 in ah three cases. This age estimate from 
population modeling of the stellar continuum is in good 
agreement with the detection of Wolf-Rayet populations, 
even though they were constrained independently. 

In Table [5] we also list the average star formation rates 
(SFR) in the past < 0.1 Gyr estimated from the stellar 
masses contained in the young stellar population divided 
by its estimated age. The inferred median SFR in our 
sample is ~ 44 Mq yr _1 (uncorrected for extinction). 
This is smaller than the median SFR ( ~ 87 Mq yr"*) 
inferre d from the IR luminosities of the Zakamsk a et al.l 
(2008) sample of 12 SDSS type 2 quasars (with a me- 
dian i[om] of 10 9 ' 4 L©, comparable to that of our sam- 
ple) assuming the calib ration of SFR-L fir for starburst 
from lKennicutH (|1998h. In pa rticular. SDSSJ0056+0032 
is also in the IZakamska et all (|2008h sample and our es- 
timate of its SFR from population synthesis (~ 36 Mq 
yr _1 ) is smaller than that inferred from its IR luminos- 
ity (~ 137 Mq yr _1 ). The discrepancy in its two SFR 
estimates suggests an extinction level of ~ 1.5 mag. 

Several caveats must be considered here. First, the 
results are all based on spectra uncorrected for internal 
reddening or extinction so that the inferred contribution 
from young stellar populations and the stellar luminosity 
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Fig. 8. — Stellar population analysis results for objects without detectable scattered quasar light. Two examples are shown here, for 
objects having the smallest (SDSSJ0056+0032) and largest (SDSSJ0950+0511) <r* measurements in our sample. The normalized spectra 
(smoothed with a 6-pixel boxcar) are shown as thin curves, and best-fit continuum models are plotted as thick gray curves. Residuals from 
the fits are displayed on the bottom along with 1-<t error spectra (dotted) and null residuals (dashed). The two fit stellar populations are 
also sh own, labelled with stellar age. Young populations (< 0.1 Gyr) contribute a considerable fraction of the total stellar luminosity. See 
3 1.3. 2 1 for more information. 
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Fig. 9. — As in Figure [8] but for an object with detected scattered quasar light as apparent from the detection of a broad component to 
H/3. This object has the highest scattered-light contribution to the monochromatic luminosity £5100 (~55%; Table [3} in our sample. The 
notation is the same as in Figure [8] the scattered light model is plotted as a long-dashed curve. See £14.3.21 for more information. 



and mass estimates could be underestimated. If scatter- 
ing is by dust, then the SED of the scattered light may be 
bluer (or redder) than the intrinsic a v = —0.44, leading 
to an overestimate (or underestimate) of the young pop- 
ulation contribution. The effects of dust in the narrow- 
line region are discussed in 95.31 Third, while we have 
subtracted scattered quasar light based on the detection 
of broad H/3 f §4.2p . there still might be some underesti- 
mated scattered light which could contaminate the stellar 
continuum if the broad component had width ^> 8000 
km s _1 . Nevertheless, the effect of this possibility on 
our results is likely to be small, considering broad-line 
quasars with such extreme line widths are quite rare. Fi- 
nally, since our targets were selected from luminous type 
2 quasars, the results on host galaxies may not necessar- 
ily apply to un-reddened luminous type 1 quasars, since 
they may trace different evolutionary phases (discussed 
below). 



4.4. Double Cores, Companions, and Extended 
Emission Line Regions 

In Figure [10] we show spatial profiles of the nine tar- 
gets constructed by collapsing the spectra over the whole 
observed wavelength range, centered on the most lumi- 
nous peak along the spatial direction. The spatial pro- 
file of the photometric standard star is also shown for 
comparison. While in many objects the central spatial 
profiles are consistent with a single point-spread function 
(PSF), those of at least two objects (SDSSJ0823+3231 
and SDSSJ0950+0511) deviate significantly from a sin- 
gle PSF and exhibit evidence for two components within 
the central 5 arcsec (what we call "double cores"). Ar- 
rows mark all the cases of double cores, companions, and 
extended emission line regions (we refer to the latter two 
cases as "companions" in general) that are physically as- 
sociated with the central components. The physical as- 
sociation is verified by determining the redshifts of com- 
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Fig. 10. — Spatial profiles of all objects in our sample and of 
the standard star (top panel marked "PSF") as measured from our 
long-slit observations. These spatial profiles have been normal- 
ized and centered on the central most luminous peaks, combined 
over the whole observed wavelength range in a single exposure for 
each object. The positions of double cores, close companions, and 
extended emission-line regions are indicated by arrows. 1 arcsec 
corresponds to 6 kpc at redshift of 0.5 in the assumed cosmology. 
The component at ~ 4 arcsec from SDSSJ0056+0032 and that 
at ~ — 7 arcsec from SDSSJ0319— 0058 have line-of-sight velocity 
offsets > 10 3 km/s relative to the central component; we do not 
consider them to be physically related companions. The spectrum 
of the neighboring component in SDSSJ0210— 1001 is too noisy to 
allow a redshift to be measured. Our long-slit observations show 
that at least four of the nine targets have double cores and close 
companions or extended emission-line regions. The profiles are 
con sist ent among different exposures for the same object. Fig- 
ure [Tl] displays 2D spectra for the two objects with double cores 
(SDSSJ0823+3231 and SDSSJ0950+0511). See Table |6] for mea- 
surements of the neighboring components. 

panions using either emission or absorption lines. The 2D 
spectra of the two objects with double cores centered on 
H/3 and [O m] along the wavelength direction are shown 
in Figure [TT] 

There are several possibilities concerning the nature of 
the companions, which include: (a) a cloud of gas which 
either is photoionized by the central nucleus and/or 
heated by shocks, or whose dust is reflecting nuclear 
light; (b) a region of active star formation; and (c) a com- 
panion galaxy which cither is merging with the central 
nucleus or is a leftover from a past interaction event, or a 




4850 4900 4950 5000 4850 4900 4950 5000 

Rest Wavelength (A) Rest Wavelength (A) 



Fig. 11. — 2D spectra of the 2 objects in our sample containing 
double cores. The spectra shown are before flux calibration and 
are uncorrected for telluric absorption, which causes the trough 
at ~ 4990 A in SDSSJ0950+0511. Contours are plotted in log- 
linear scales; the lowest and highest levels are 

10 i.5*0 ( 10 1.25) and 

10 4.50 ( 10 3-50) counts witn 

a spacing of 0.50 (0.25) dex, amounting 
to 7 (10) levels in total for SDSSJ0823+3231 (SDSSJ0950+0511). 
Spectra have been centered on the most luminous peaks along the 
spatial direction. The position of the second central component 
is marked with a red dashed line. Dotted lines mark the position 
of the aperture used to extract ID spectra. In SDSSJ0823+3231, 
the second central component is most prominent in emission lines 
whereas in SDSSJ0950+0511 it is most prominent in the contin- 
uum. There are additional companions and extended emission line 
regions with spatial offsets ASlos^ i5 arcsec. The LOS veloc- 
ity offsets of these components are indicated on the plot. The 
two neighbor components of SDSSJ0950+0511 have clear continua 
and much smaller [Om]/H/3 ratios (> 1) than that in the cen- 
tral nucleus; these ratios are more characteristic of star-forming H 
II regions. In contrast, the two components of SDSSJ0823+3231 
exhibit almost the same [Our]/H/3, [O iii]/[0 n], and [Ne v]/H/3 ra- 
tios as those of the nucleus, characteristic of AGN excitation. The 
component at a projected spatial offset of ~ — 5 arcsec has a clear 
continuum detection (although it is not apparent in the plot due 
to the scaling) whereas that at ~ 5 arcsec has a weak continuum 
barely detected redward of ~ 5100 A. 

small satellite which is not interacting with the nucleus at 
all (at the time of observations). Discriminating among 
these possibilities is useful to assess the significance of 
galaxy interactions in triggering luminous quasar activ- 
ity. We list the projected spatial offsets, line-of-sight ve- 
locity offsets, SDSS magnitudes, and several diagnostic 
emission-line measurements of the companions and dou- 
ble cores in Table El A companion to SDSS J0056+0032 
and a companion to SDSSJ0319— 0058 with relative ra- 
dial velocities > 1000 km/s are not listed. 

The companions in SDSSJ0319-0058 and in 
SDSSJ0950+0511 are likely to be regions of active 
star formation, given their diagnostic e mission line 
ratio s [On]/[Om] and/or H/3/[Om] (e.g., lOsterbrockl 
119891 ). For the companion in SDSSJ0319-0058, there 
is no evidence for disturbance or interaction with the 
central nucleus in its 2-D spectrum (e.g., multiple 
components at different velocities or tidal tails). In the 
double-core system SDSSJ0950+0511, it is likely (see 
below) that galaxy interaction is responsible for both 
the quasar and the young starburst activity, and that 
the luminous type 2 quasar is at an early stage of inter- 
action. The companion at ~ —32 kpc has a line-of-sight 
(LOS) velocity offset (relative to the more luminous 
core) similar to that of the less luminous core (Figure 
[PT1 right panel). This suggests that the companion is 
physically associated with the less luminous core and 
might have been stripped off as a result of the merger in 
the nucleus (which is marginally resolved in this case) . 

The companion in SDSSJ0943+3456 is an absorption- 
line galaxy. It is likely to be a relic core from a past 
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TABLE 6 

Companions, Double Cores and Extended Emission-line Regions 



Target Name . (1) SDSSJ0319 SDSSJ0823 SDSSJ0943 SDSSJ0950 



As (kpc) (2) 


-18 


5 -25 


30 


46 


-9 -32 


37 


A„ (km/s) (3) 


-198 


-124 


-170 


-96 


-383 -452 


87 


r (4) 




22.75 




21.30 


22.14 




EW [0ln] (A) (5) 


4.5 


682 


52.5 




15.9 


11.2 


FWHM [0 ni] (km/s) (6) 


167 


212 


243 




332 


640 


log(i[o m] A©) ■•■ (7) 


7.45 


9.08 


7.88 




7.93 


7.43 


[On]/[Om] (8) 




0.10 


0.43 




1.9 


1.7 


H/3/[0 in] (9) 


1.8 


0.12 


0.09 




0.48 


0.15 



Note. — 
Row(l): Abbreviated target name. 

Row(2): Projected spatial offset (physical distance) with respect to the central nucleus. The plus or minus sign follows the definition of Figurc flQI 
Row(3): Linc-of-sight velocity shift relative to the central nucleus. 

Row(4): r-band model magnitude from SDSS when available. Several components do not have SDSS photometry because they are too faint. 

Row(5): Rest-frame equivalent widths of [O ill] A5007 measured for the components containing line emission. The spectrum of the component in SDSSJ0943 

is free of emission lines. The double cores of SDSSJ0823 and SDSSJ0950 are not well spatially resolved (namely within the central 5 arcsec). 

Row(6): Full widths at half maximum of [O ill] A5007 corrected for instrumental resolution. 

Row(7): Line fluxes of [O in] A5007. 

Row(S): Emission-line ratio — [° n ] AA 3727,3729 c } larac ^ er ^ zm o- ^ nc ionization parameter. There is no measurement for the component in SDSSJ0319 since 

v ' F [OIII]A5007 
its [O n] A3727 line overlaps the CMOS chip gap. 

Row(9): Emission-line ratio -p* — . 

*[OIII]A5007 



interaction event with the central galaxy, because it has 
a best-fit stellar age of 5 Gyr from stellar population 
synthesis, which is the same as the older population in 
the circum- nucleus stellar component (Table[5]). We have 
detected several stellar absorption features (including Ca 
K&H, the G band, Mgife, and Fe 5270) in its spectrum 
and measured its <r* to be 295 ± 9 km/s. 

According to the diagnostic line ratios [On]/[Om] and 
H/3/[0 m] the companions in SDSSJ0823+3231 are likely 
to be gas clouds either photoionized by the central nu- 
cleus and/or heated by shocks, or their dust is reflecting 
quasar light, rather than regions of active star formation. 
They are unlikely to be galaxies, as there are no stellar 
features detected in the continua. The non-stellar con- 
tinua are most likely scattered quasar light. The emis- 
sion lines seem to be from gas photoionized by the central 
nucleus. The absence of broad lines, and the large EWs 
of the narrow emission lines suggest that these lines are 
not due to scattering from gas clouds within the ion- 
izing cones; scattered light from gas clouds outside the 
ionizing cones are unlikely to produce the bulk of the 
observed luminosities of the companions. Shock heating 
also seems unlikely as the FWHMs of the emission lines 
are small. Since SDSSJ0823+3231 is a double-core sys- 
tem, these gas clouds might be gas shreds resulting from 
the on-going merger. 

In summary, our long-slit spectroscopy reveals that at 
least four of our nine targets contain double cores and/or 
physically associated companions. Three of the four ob- 
jects show evidence for galaxy interactions which may 
be responsible for the quasar and starburst activity. The 
companions in these systems have various origins, includ- 
ing a leftover galaxy core dominated by old stellar popu- 
lations, merging galaxies with active star formation, and 
perhaps gas clouds shredded by the merger. While we 
oriented the slit to cover potential companions seen in 
the imaging data, our spectra do not resolve scales < 1 
arcsec so that the interaction fraction we find is a lower 
limit to the true value. 



5. IMPLICATIONS AND DISCUSSION 

In this section we present implications and discussion 
on our results. We infer black hole masses and Eddington 
ratios in §5.1[ and discuss the origins of the blue continua 
in luminous type 2 quasars in A5.2I with a prescription 
for the scattered quasar light in i )5.2.1l and the star burst- 
quasar link in §5.2.21 Extinction and its effects on our 
results are provided in t )5.3l In t )5.4[ we compare our 
results with those from other quasar host-galaxy studies, 
both obscured ( §5.4. 1|) and unobscured ( §5.4.211 . 

5.1. Black Hole Mass and Eddington Ratio 

We infer black hole masses from the measured stel- 
lar velocity dispersions ( §4.ip based on the correlation 
between bulge stellar velocity dispersion and dynami- 
cal black hole masses obse r ved in local inactive g alaxies 
(|Ferrarese k Merrittl 120001: iGebhardt et all [2000) . using 
the calibration of iTremaine et al.l (|2002f ) The calibra- 
tion for local inactive galaxies do es not necessarily d i- 
rectly apply to active galaxies (e.g. [Greene XHc][2006bl) . 
considering that black hole s in active galaxies are still in 
growth (e.g- lHo et al.l[2 008). Further more, while our tar- 
gets only have moderate redshifts z ~ 0.5, there could be 
non- negligible redshift evolution in the Mb h-c* relation 
(e.g. IWoo et"aLll2008l : ICanalizo et al.ll2008l ). For objects 
with detected broad H/3, we examine this redshift evo- 
lution in S I5.1.2I We estimate accretion rates bearing in 
mind these uncertainties. Adopting £[om] as a proxy 
for quasar intrinsic power (see E I5.1.ip . we obtain quasar 
bolometric luminosities and Eddington ratios. 

Black hole masses, quasar bolometric luminosities, and 
the inferred Eddington ratios are listed in Table [2] and 
displayed in Figure fT2l Our sample has a median black 
hole mass of 1O 8 8 M0 and a median Eddington ratio 
of - 0.7. As cautioned in gTj SDSSJ0823+3231 and 
SDSSJ0950+0511 have resolved double cores in their cen- 
tral 5 arcsec which could bias cr* towards larger values, 
resulting in overestimated black hole masses and under- 
estimated Eddington ratios for the central nuclei. Tak- 
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Fig. 12. — Quasar bolometric luminosity -LboIi black hole mass 
Mbh, and Eddington ratio iflol/^Edd- Left: LboI versus Mbh- 
The loci with ^Bol/^Edd equal to 10, 1, and 0.1 are plotted as dot- 
ted, solid and dashed lines. The two objects indicated with open 
circles h ave marginally resolved double cores within the central 5 
arcsec ( 34.4)1 and the measurements represent luminosity- weighted 
values of both components. Right: Distribution of Eddington ra- 
tios. Our sample has a medianLBoi/^Edd of 0.7. We list the 
measurements in Table [2] See 35.1l for more discussion. 

ing into account all the associated uncertainties, almost 
all of our targets are accreting at higher than 10% of 
the Eddington rates. The two targets estimated to have 
super-Eddington ratios both have Wolf-Rayet detections. 
The Eddington ratios we fi nd overlap with th ose of SDSS 
type 1 quasars studied by IS hen et al.l (|2008f ) with com- 
par able bolometric lum inosities at similar redshifts (also 
sec iGreene et al1l200l . 



The contamination on L 



[Oiii] 



from star formation 



in our targets seems negligible given their [On] /[Om] 
emission-line flux ratios listed in Table [7] (e.g., 
iFerland fc Osterbrockl 1986; Hd l2005l) . This statement 
still holds in the two super-Eddington targets. 

5.1.1. Bolometric Correction for Type 2 Quasars 

We estimate bolometric luminosities LboI of type 2 
quasars based on L\o UI i. First, we infer the intrinsic 
-^2500 from L[om] using the calibration by iReves et al.l 
(2008) of the M 25 oo-£[Oiii] relation observed in SDSS 
type 1 quasars. Then we estimate LboI from M2500 
using the bolometr i c correction for ty pe 1 quasars 
(|Marconi et al.ll2004t iRichards et al.ll2006h and estimate 
the systematic uncertainty of the bolometric correction 
(the M2500 to LboI conversion) to be 0.1-0.2 dex in the 
relevant luminosity range (see below). The final calibra- 
tion obtained is given by 



log 



^Bol 

La, 



0.99 x log 



J [0 in] 



+ 3.5, (1) 



with a total l-a uncertainty of 0.5 dex on log(LBoi/-^o) 
at LboI ~ 10 13 L Q . The total uncertainty of LboI con- 
volves the errors propagated from the measurement error 
of I/[ Q m |, the 0.3 6 dex scatter of the M25oo-£[om] rela- 
tion (IReves et al. 2008), the 0.05 dex scatter on Lb/LboI 
(|Marconi et al.ll2004h . and the additional 0.1-0.2 dex sys- 
tematic uncertainty (depending on £[Om]) from bolo- 
metric corrections using different estimates. The total 
uncertainty is dominated by the intrinsic scatter of the 
M 2 5oo-i[Om] relation. 

The systematic uncertainty of the bolometric correc- 
tion for type 1 quasars is estimated by comparing three 
correction approaches: (1) converting from M2Rn n to Lb 
assuming a v = -0.44 (|Vanden Berk et alj[2001j) and to 
^Boi with the luminosity-dependent £?-band bolometric 



correction of iMarconi et al.l (|2004D ; (2) converting from 
M2500 to LboI usi ng the luminosi t y-inde pendent bolo- 
metric correction of lRichards et a l. (2006); and (3) same 
as (2), but in the calculation of the bolometric correc- 
tion, we extrapolate the 5000-10000 A quasar spectrum 
into the mid-IR to avoid doub l e counting of re-ra diated 
emission (jMarconi et all l2004t IReves et all I2008D . At 
L[owi\ > 10 9 '°£q, method (2) gives results <0.2 dex 
higher than those of method (1), while method (3) gives 
results <0.2 dex lower than those of method (1). We 
take the results from method (1) as our baseline values 
and the differences among the three as our estimate of 
the systematic uncertainties. 

5.1.2. Comparing black hole Mass Estimates Based on 
a* and Scattered Quasar Light 

For the subset of targets which have broad-H/3 detec- 
tion ( §4.2p . we could estimate virial black hole masses 
using the FWHMs from broad H/3 measurements and 
the restored intrinsic quasar continuum luminosity £5100 
(Table [3]) converted from L\q m | using the M 2 5oo-£[Oiii] 



relatio n of IReves et al.l ( 20081 and assuming a v = 
-0.44 ()Vanden Berk et alJ I^OOl). Although it depends 
on the uncertain H/3 FWHM mass inferences (e.g., 
IMarconi et al.l 120091 ). and on the bolometric corrections 
for type 2 quasars, which again have a large systematic 
uncertainty, the scattered broad line approach allows us 
to independently calibrate the estimation of black hole 
mass from host galaxy properties (bulge stellar velocity 
dispersions and luminosities) in type 2 quasars. Results 
on the black hole m ass -Mgjj using the virial formula (e.g. 
IGreene fc Hcl2005l) are listed in Tabled The uncertainty 
of is dominated by that propagated from the quasar 
intrinsic power. While the two approaches for estimat- 
ing black hole mass are based on completely different 
assumptions and have separate systematic uncertainties, 
the virial estimates based on scattered quasar light are 
consistent with cr* -based black hole masses within the es- 
timated uncertainties. This general agreement lends fur- 
ther support to the robustness of both our scattered-light 
and a* measurements, at least for the small subset of tar- 
gets having broad-H/3 detections. Based on our current 
data, we cannot draw any firm conclusion on redshift evo- 
lution considering the large uncertainties and the small 
sample size. While the associated uncertainties are large 
in practice, the approach offers the possibility to test the 



at low redshifts (jLauer et al 



ge where 
U l200l . 



5.2. Origins of the Blue Continuum in Luminous Type 
2 Quasars 



We find that quasar scattered light fi j4.2[) and starlight 
from young populations f ^4.3.2p both contribute consid- 
erably to the blue continuum observed in luminous type 
2 quasars. Their relative importance varies among differ- 
ent objects, with scattered light dominating in some cases 
(e.g., SDSSJ0823+3231), and massive stars dominating 
in others (e.g., SDSSJ0056+0032). In this section, we 
first present the dependence of scattered light on L[ 0m ], 
and then discuss the implications of the prevalence of 
young stellar populations in luminous type 2 quasars. 

5.2.1. A Prescription for Scattered Light 
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The scattered light contamination in low luminosity 
type 2 AGN rarely exceeds 5% in the rest-frame opti- 
cal (e.g.lCid Fernandes fc Terlevichlll995t ISchmitt et all 
119991 Kauff mann et alll2003fh ~ However, we find that it 
can be signi ficant when L[Q m ] is large ( §4.21) . confirming 
the results of lZakamska et all (j2005Ll2006f ) based on spec- 
tropolarimetry and HST imaging that a significant frac- 
tion of continuum emission can be due to scattered light 
in luminous type 2 quasars. The ratios L|^Q t tered /L§!j'g 
for our nine targets range from < 5% to 55% (Table [SJ; 
combined with the estimated quasar intrinsic luminosity 
from L[q m ] ( §5.1.1|) . the inferred scattered light fractions 

L tToQ CIcd / L ?im for our ninc targets range from < 0.2% 
to 0.4%. In view of the high incidence of scattered light 
in luminous quasars, its quantification is particularly im- 
portant for robust determination of stellar populations 



Figure [TBI displays both scattered light luminosity and 
its fractional contribution as a function of Lroml- Also 
shown are measurements compiled from the literature for 
comparison (see figure caption for details). The best-fit 
linear models to all the data shown are given by 
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While there is substantial scatter in the above relations 
and the amount of scattered light received could vary due 
to different coverage of the scattering cone by different 
observations even at a given observation angle, the em- 
pirical models are useful for a quick estimate of scattered 
light given any L[om]- We caution that the scattered 
light fraction estimated here (within the slit) is a lower 
bound to the total scattered light fraction, since there 
could be considerable amount of scattered light outside 
our slit. 

5.2.2. Young Starburst 

After the removal of scattered light, there is still con- 
siderable blue light in the spectra of our targets even 
before reddening correction. We found significant con- 
tributions (> 50% of scattered-light subtracted £5100) 
from very young populations (< 25 Myr) in four out of 
nine of our objects; eight of our nine targets have > 20% 
of L5100 contributions from very young stellar popula- 
tions (Table [5]). In six of these objects very young stars 
account for < 5% of the stellar mass. 

Considering the short duration of the WR phase (2-5 
Myr after a burst of star formation), the fact that > 30% 
of our objects show WR features suggests that if the star- 
burst and Eddington-rate feeding to the SMBH proceed 
coevally, the Eddington-rate quasar activity directly as- 
sociated with the starburst lasts < 10 Myr, assuming 
type 2s do not evolve into type Is. While a more careful 
estimate for the timescale can be obtained by assum- 
ing models for the quasar light curve and the evolving 
circum-nuclear obscuration, the crude estimate here is 
in agree ment with that based on black hole demograph- 
ics (e.g. lYu fc Tremaindl2002l ). Enhancing the statistical 
significance of our results will require studying a larger 
sample of Eddington-rate type 2 quasars. The general 
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Fig. 13. — Scattered light contamination as a function of £[om] 
in the on-nucleus spectra of luminous type 2 quasars. For the 
nine targets in this work the estimates are based on broad-H/3 
detections ( 94.20 and the correlation between H/3 and continuum 
luminosities observed in unobscured quasars (Greene & Hol l2005T ). 
Also displayed in comparison are scattered-light estimates com- 
piled from the literature in cluding: (1) IRAS 09104+4109, based 
on spectropolarimetry from Tran et al. (2000), assuming an intrin- 
sic continuum polarization of 20%; (2) 12 SDSS type 2 quasars, 
based on continuum decomposition and broad-band (downward 
triangles for yellow band and upward triangles for g reen band) 
polarization measurements from Zakamska et al. (2005), assuming 
an intrinsic continuum polarization of 20%; and (3) upper limits 
on scattered light for 15 SDSS type 2 quasars based on broad- 
H/3 measurements (Greene et al. 2009;). One of our Gemini tar- 
gets (SDSSJ0157— 053) has also been independently studied by 
I Greene et al. (2009) and its two measurements are connected by 
a straight line. The long dashed line presents the best-fit linear 
model to all the measurements shown, and the dashed lines denote 
1-<T scatter. The models are given in Equation [2] 

approach of using WR features as a clock demonstrated 
with our pilot sample here is another way to observation- 
ally estimate the Eddington-rate quasar duty cycle. 

Our data do not have the spatial resolution to de- 
termine if these starburst events are in the immediate 
circum-nuclear regions or are occurring on much larger 
galaxy scales. Nevertheless, it appears that galaxy inter- 
actions are responsible for or coincident with both the 
quasar and the starburst activity, at least in the double- 
core systems fi j4.4[) . 

5.3. Extinction and Its Effects on Our Results 

Most of the results above are based on spectra un- 
corrected for dust extinction in the galaxies themselves. 
In Table [7] we list reddening estimates in the narrow- 
line regions from the intensity ratios tPy/H/? and H5/H/3 
using the Balmer decrement method, assuming the in- 
trinsic case B values of 0.466 a nd 0.256 respectively for 
T = 10 4 K (|Osterbrockill989D . the extinction curve of 
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TABLE 7 

Emission-line Ratios and Color Excess 



Target Name fgj B(B - V) E(B-V) Q 

(1) (2) (3) (4) (5) (6) 



SDSSJ0056+0032 0.352 ± 0.004 0.142 ±0.003 0.54 ± 0.02 0.77 ±0.03 0.14 ±0.01 

SDSSJ0134±0014 0.400 ± 0.005 0.203 ± 0.004 0.29 ± 0.02 0.30 ± 0.03 0.09 ± 0.01 

SDSSJ0157-0053 0.455 ± 0.004 0.245 ± 0.003 0.05 ± 0.02 0.06 ± 0.02 0.15 ±0.01 

SDSSJ0210-1001 0.659 ±0.005 0.342 ± 0.003 -0.67 ±0.01 -0.38 ± 0.01 0.17 ±0.01 

SDSSJ0319-0058 0.563 ± 0.006 0.277 ± 0.004 -0.36 ± 0.02 -0.10 ±0.02 0.14 ±0.05 

SDSSJ0801±4412 0.469 ± 0.008 0.242 ± 0.005 -0.01 ± 0.03 0.07 ± 0.03 0.13 ±0.01 

SDSSJ0823±3231 0.400 ± 0.004 0.224 ± 0.003 0.29 ± 0.02 0.18 ±0.02 0.18 ±0.01 

SDSSJ0943±3456 0.472 ± 0.005 0.199 ±0.003 -0.02 ± 0.02 0.33 ± 0.02 0.10 ±0.01 

SDSSJ0950±0511 0.383 ± 0.007 0.104 ±0.005 0.38 ± 0.04 1.18 ±0.06 0.18 ±0.01 



Note. — 

Col. (2), (3): Emission-line ratios measured over continuum-subtracted GMOS spectra uncorrected for reddening. The quoted 
1-a uncertainties were derived based on error spectra. 

Col. (4). (5): Color excess estimated from the intensity ratios (for Col. 4) and j?^* 5 (for Col. 5), assuming the intrinsic 

case B values of 0.466 and 0.256 respectively for T — 10 4 K ( Ostcrbrock 1989) and the extinction curve o f Cardclli ct al. (1989) 
with i?v — 3.1. For most of our targets (7 out of 9), the two reddening estimates do not agree with each other and some are 
negative, indicating deviations from standard reddening laws and the assumed foreground dust-screen model. 

Col. (6): Emission-line ratio ^ °f ^ \ \ qq^ 2 9 characterizing the ionization parameter. The [O n] A3727 line of SDSSJ0319— 0058 

overlaps the GMOS chip gap and its listed value is measured from SDSS spectra. 



ICardelli et alj (|1989f ) with Ry = 3.1, and a model of 
a foreground obsc uring screen. However, as found by 
IReves et all (|2008ft in a subset of the parent quasar sam- 
ple having -Lroin] > 10 90 L@, the intensity ratios Ha/H/3 
and H7/H/3 do not obey standard reddening laws, nor 
can they be described by a simple dust screen. We find 
a similar situation in our sample by comparing the in- 
tensity ratios H7/H/3 and H5/H/3, seven out of nine of 
which have significantly different reddening estimates for 
the same object (Table [7]). 

Despite this caveat, the estimated color excesses E(B — 
V) for our 9 targets range from —0.67 (—0.38) to 0.54 
(1.18) with a median of 0.05 (0.18) according to H7/H/? 
(H<5/H/3). The color excess estimated for the ULIRG, 
IRAS 09104±4109, bv lTran et aD (j200l is E(B - V) = 
0.24, also using the Balmer decrement combining the 
intensity ratios Hy/H/3 and H(5/H/3, a value very simi- 
lar to those found in our sample. However, the small 
color excesses derived using the Balmer decrement do 
not necessarily represent the true dust content of these 
objects; IRAS 09104+ 4109 is hyper- luminous in the in- 
frared (Lir, > 10 13 i©; Klcinma nn et al.lll988h yet has a 
small color excess according to the Balmer decrement. 
This just suggests that most of the dust is concentrated 
on scales smaller than the Balmer lines are emitted. High 
extinction with little reddening could arise in objects 
with very patchy and optically-thick dust clouds. 

In view of the uncertainties in extinction geometry 
(e.g., §6.2 of IReves et all [20081. we do not apply cor- 
rections but rather keep in mind the potential effects on 
our results. The unaccounted extinction would also bias 
emission-line luminosity measurements smaller, which 
means the quasar luminosities and Eddington ratios 
ft j5.ll Figure [T2"|) could be underestimated. The unac- 
counted reddening would bias the stellar ages higher and 
the estimates of the contribution from young populations 
lower (§1X21). 



5.4.1. Host Galaxy Studies of Type 2s 

We now compare our results with other host galaxy 
studies of type 2 AGN. The present work is the first 
deep optical spectroscopic study of luminous type 2 
quasars with Lrom] > 10 9 ' 3 L Q . A number of previ- 
ous groups have found that the presence of young stel- 
lar populations is a general prop erty of type 2 AGN 
with high [O m] luminosities ( e.g. [S chmitt et al. 119991 : 
i Gonzalez Deleado et al.l I2001L IKauffmann et alJ 120031 : 
ICid Fernandes et alj 12004 IWild et al.ll2007ft . However, 
the typical ages of young stellar populations (<0.1 Gyr) 
we find are much sm aller than the pos t -starb urst ages 
(1-2 Gyr) seen in the [Kauffman n et all (|2003l ) subsam- 
ple with Lp nr i/iQ ~ 10 7 -°-10 8 - 5 (reddening corrected). 
This result suggests that the luminosity dependence of 
host galaxy stellar age a nd post-starburst fraction seen 
at lower luminosities (e .g. lHeckmanlll980l : iHo et al.ll2003t 
IKauffmann et all l2003f ) extends to the most luminous 
quasars. Within our luminous quasar sample, we do not 
see a dependence of stellar age on iromlj but the scatter 
of the correlation (if any) should be large, and our sam- 
ple may simply be too small to see such a correlation. 
A much larger sample is needed to study the luminosity 
dependence of host galaxy properties in luminous type 2 
quasars. 

Most intriguingly, we detect WR populations in a third 
of our sample; WR feat ures are much less often seen a t 
lower luminosities fe.g. IGonzalez Delgado et alj 12001). 
The fraction should be treated as a lower limit since 
it does not account for regions so obscured that they 
do not contribute to the optical emiss ion. If starburst 
and q uasar activity peak coevally (e.g. iDi Matteo et al.1 
l2005h . the chance of detecting the WR phase would be 
much larger in luminous quasars close to the peak activ- 
ity than in less luminous quasars which are in their later 
stages of evolution, after WR features fade away and 
intermediate-age post-starburst signatures take over. 



5.4. Comparison with Other Quasar Host Galaxy 
Studies 



5.4.2. Host Galaxy Studies of Type Is 



18 



LIU ET AL. 



We now compare our results with those from host 
galaxy studies of type 1 quasars and address the issue 
of whether t ype 2s and type Is follow a temporal se - 
quence (e.g.. ISanders et al.l IT988; Ho pkins et aLll2006f ). 
Our type 2 targets have quasar intrinsic luminosities of 
My < —26 mag ( §5.1.1|) . There is virtually no spectro- 
scopic information on host gal axies of type 1 quasar s at 
this luminosity scale (but see IBoroson &; OkdTl982l . for 
3 C 48). 

iFlovd et all (l200l imaged 17 type 1 quasars with 
-28 < My < -24 mag and z ~ 0.4, and found them 
to be massive bulge-dominated galaxies. The Edding- 
ton ratios these authors find are comparable to those of 
our targ ets, despite the very different systematic uncer- 
tainties. lUrrutia et all (|2008[ ) studied 13 reddened type 

1 quasars with —26 < Mg < —24 mag and z ~ 0.7, 
and found that 11 of them show evidence for recent in- 
teractions, a much high er fraction than tha t found by 
IFlovd et alJ (j20?)l and iGuvon et alJ (pOfil . While it 
is possible that deeper high-resolution im aging will re- 
veal w eak interaction features (see below). lUrrutia et al.l 
(2008) suggest that this inconsistency can be explained 
in the ev olutionary scenario fo r quasar obscuration pro- 
posed by ISanders et al.l (|1988j ). in which dust-reddened 
type 1 quasars trace an earlier evolutionary phase than 
ordinary type 1 quasars, and thereby are more prone to 
signs of a merger. 

If the evolutionary argument is true, we should also see 
a hig h interaction fra c tion i n our obscured quasar sam- 
ple. IZakamsk a et alJ (|2006l ) found almost half in their 
obscured quasar sample show signs of mergers or inter- 
actions. Only two of our targets show double cores (S j4.4[) , 
but the stellar population results clearly suggest that our 
targets trace an early stage of galaxy interactions if they 
are responsible for the young starburst activity. 

While the strict AGN unification model states that 
type 2s are the same as type Is other than the differ- 
ence in our viewing angle (e.g., lAntonuccll ll~993) . it is 
likely that orientation and evolutionary effects are both 
at work, and both need to be considered to interpret the 
selection of a type 2 quasar sample. If the opening angle 
increases as an object proceeds from the reddened to the 
ordinary quasar phase and the ordinary phase does not 
last much longer than the reddened phase (a statement 
that needs to be tested), both of which have ~ Eddington 
accretion rates, then by construction (averaged over both 
orientation and time), the majority in a sample of lumi- 
nous type 2 quasars would be the obscured counterpart of 
reddened luminous type 1 quasars, with only a minority 
resembling regulars. In at least some type 2 quasars, ex- 
tinction processes are observed to operate both at small 
scales due to circum-nuclear d usty tori, and at galaxy- 
scale s due to star formation (e.g.[ Martmcz -Sansigre et al.l 
120061: iRlgbv et alj|200l iLacv et al.ll2007l) . Using Spitzer 
spectra. IZakamsk a et alJ ( |2008l) show that SDSS type 

2 quasars have the highest star-formation rate among 
all quasar samples except for other samples of type 2 
quasars. Luminous type 2 quasars, the majority of which 
resemble reddened type 1 quasars, should be sensitive to 
the early stages of quasar evolution, and hence are ideal 
for quantifying the role played by galaxy interactions. 

Quasars with lower luminosities (—25.5 < My < 
—23.5 mag) and with smaller Eddington ratios (0.01 < 
L/LEdd < 0.1) almost all reside in massive ellipticals 



with relaxed light profiles (|Dunlop et al.ll2003[). Deeper 
imagin g studies of a few objects in the iDunlop et al.l 
(2003) sample have revealed shells and ripples which are 
interpreted as merger relics, with ages of several hun- 
dred of Myr to 1 Gyr estimated based on merger simu- 
lations (jCanalizo et alJ 120071 iBennert et all 12008). con- 
sistent wit h spectroscopically d etermined post-starburst 
ages (e.g. ICanalizo et all 120071 1.4 Gyr). While shells 
and ripples are telltale signatures of galaxy interactions, 
these systems are observed at a late stage of mergers, 
similar to th e small subset associat ed with mergers in 
type 2 AGN (jKauffmann et al J 120031). The typical post- 
starburst ages of 1 -2 Gyr (see also Uahnke et alJ l2007t 
iLetawe et al.ll2007h are significantly larger than those of 
our targets (< 0.1 Gyr). 

In summary, our results combined with those of other 
quasar host galaxies spanning different ranges of quasar 
luminosity and Eddington ratio may indicate that host 
galaxy stellar population properties depend on both 
quasar luminosity (or more fundamentally Eddington 
ratio) and obscuration. This is because the two fac- 
tors both characterize the phase of quasar evolution and 
by extension that of stellar evolution if the growth of 
SMBHs and stellar bulges are coupled. Clearly more 
spectroscopic data and images with higher resolution are 
needed to better understand the host galaxies of both 
type 1 and type 2 quasars with luminosities comparable 
to luminous z > 2 quasars. It is important to quantify 
the fraction of objects containing WR populations in lu- 
minous type 1 quasars and to compare with our results 
for luminous type 2 quasars, although doing so will be 
quite challenging because of the brightness of the quasar 
continuum. 

6. SUMMARY 

We present deep Gemini GMOS optical spectroscopy of 
nine luminous type 2 quasars at redshift s z ~ 0.5, drawn 
from the S PSS type 2 quasar sample of IZakamska et all 
(|2003h and lReves etail (|2008l ). Our targets were selected 
to have high intrinsic luminosities (My < —26; £bo1 ~ 
1O 13 L0) inferred by the [Om] luminosity. 

Our main findings are as follows: 

1. We unambiguously detect WR populations in three 
of our nine objects. These very young post- 
starburst populations are independently confirmed 
by population modeling of the stellar continuum: 
in all three targets with WR detections, a 5-Myr 
post-starburst population contributes substantially 
(> 50% of £5100) to stellar continuum luminosities. 
In five of the six targets without WR detections, 
population modeling shows the presence of consid- 
erable (> 20% of L5100) young (< 0.1 Gyr) post- 
starburst populations. We estimate a median SFR 
of our sample of ~ 44 M© yr" 1 (uncorrected for 
extinction) . 

2. We have subtracted continuum scattered light be- 
fore performing our stellar population synthesis 
analysis. The scattered light is quantified by 
measuring a broad H/3 component underlying the 
strong narrow H/3 line. The broad H/3 component is 
detected in four targets, and the inferred scattered 
light contaminates 30-60% to the total observed 
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continuum at 5100 A. The scattering fraction 
(within the slit) is estimated to be < 0.4% at 5100 
A in all our objects, while we caution it is a lower 
bound to the total scattered light fraction, since 
there could be considerable amount of scattered 
light outside our slit. We provide a prescription for 
scattered light as a function of £[om] (Equation^ 
which can be used to get a quick but rough estimate 
of the contamination level. 

3. We obtain stellar velocity dispersion measurements 
er* for the host galaxies from the G band stellar 
absorption feature, which range from 120 to 420 
km/s in our sample. These correspond to black 
hole masses of ~ 10 7 ' 2 -10 9 4 M Q with a median 
va lue of 10 8 ' 8 M(T), a ssumi ng the Mbh-c calibration 
of iTremaine et ail (|2002l ). Combined with the es- 
timated intrinsic quasar power from L[q mi (Equa- 
tion[lJ, the inferred Eddington ratios are ^0.1-12, 
with a median of 0.7 and a scatter consistent with 
measurement uncertainties. 

4. Our long-slit spectroscopy shows that four of our 
nine targets contain double cores and/or physically 
associated companions. At least in three systems, 
a galaxy interaction might be responsible for both 
the quasar and the young starburst activity, and 
that luminous type 2 quasars with Eddington ra- 
tios close to unity perhaps trace an early stage of 
interaction. Our results combined with those from 
other quasar host galaxy studies may suggest that 
host galaxy stellar population properties depend on 
both luminosity and obscuration which character- 
ize the phase of quasar evolution. 

By targeting the obscured population, we are able to 
study stellar populations in quasars that are as luminous 
as the majority of z > 2 quasars. We are directly seeing 
the youngest stellar populations in quasar hosts, pro- 
viding strong evidence for a direct link between on-going 
starburst and luminous quasar activity. However, our pi- 



lot sample is too small to study the correlations between 
quasar and host galaxy properties, to probe the physi- 
cal links and mutual influence of starburst and luminous 
quasar activity. Correlation studies of AGN and their 
hosts at lower lu minosities have yielded extensive phys- 
ical i nsight (e.g., Kauffmann ct al. 2003; Hcckma n et al.l 
l2004h but are only enabled with the assembly of statis- 
tically large samples. We are going to carry out more 
observations in the optical to do such correlation studies 
in luminous type 2 quasars at z ~ 0.5. We look forward 
to the next generation of ground-based n ear-IR multi- 
object spectrographs (e.g.. lBeH~et al. 2009]) that will pro- 
vide statistical samples of luminous obscured quasars at 
high redshifts to better understand the coupled growths 
of stellar populations and SMBHs in the early universe. 
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